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Abstract 

We consider spherically symmetric, asymptotically flat space-times 
with a collisionless gas as matter model. Many properties of the result- 
ing Vlasov-Einstein system are not yet accessible by purely analytical 
means. We present a discretized version of this system which is suitable 
for numerical implementation and is based on the particle-in-cell tech- 
nique. Convergence of the resulting approximate solutions to the exact 
solution is proven and error bounds are provided. 

1 Introduction 

The properties of space-times filled with matter are of considerable interest 
in general relativity. The analytical and numerical feasibility as well as the 
predictions of the models depend to a large extent on the choice of matter model. 
A model for which considerable progress has been made is the collisionless gas. 
We consider a large ensemble of particles, all of which are assumed to have the 
same rest mass, normalized to unity, and to move forward in time. Therefore, 
their number density / is a non-negative function supported on the mass shell 

PM:={ 5a/s pV = -l s P°>0}, 

a sub-manifold of the tangent bundle TM of the space-time manifold M with 
metric g a p. We use coordinates (t,x a ) with zero shift and corresponding canon- 
ical momenta p a ; Greek indices always run from to 3, and Latin ones from 1 
to 3. On the mass shell PM the variable p° becomes a function of the remaining 
variables (t,x a ,p b ): 

P° = V-9 00 \/l + 9abP a P b . 
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The Vlasov-Einstein system which governs the time evolution of such a fully 
relativistic, self-gravitating collisionless gas now reads 

G al3 = 8irT al3 , 

J -Pa 
where are the Christoffel symbols, \g\ denotes the determinant of the metric, 
G a/3 the Einstein tensor, and T a<3 is the energy-momentum tensor. All physical 
constants, in particular the speed of light, are set to unity. We refer to jf] for a 
review of results on this system. Note that the characteristic system 

x a =p a /p°, p a =rf 3 y P y P 

of the Vlasov equation are precisely the geodesic equations in the given metric, 
written in coordinate time t. The characteristic flow leaves the mass shell PM 
invariant. 

Little is known about the behavior of solutions to the Vlasov-Einstein system 
in full generality. In the present investigation we restrict ourselves to the spher- 
ically symmetric and asymptotically flat case. Using Schwarzschild coordinates 
the metric is taken to be 

ds 2 = -e 2 ^dt 2 + e 2X ^dr 2 + r 2 (d9 2 + sin 2 6d<p 2 ), 

where t € 1R, r > 0, 8 €E [0, ir] , ip £ [0, 2ir) . Asymptotic flatness is then expressed as 

lim X(t,r)= lim n(t,r) = 0, tGJR, (1.1) 

r — >oc r — >oc 

and a regular centre is guaranteed by 

A(t,0) = 0, teTR. (1.2) 

To formulate the system it is convenient to let x = r(sin#cos</9,sin0sin</3,cos#) 
and, instead of using the corresponding canonical momenta to coordinatize the 
mass shell over the space-time point (t,x) use frame components 

r r 

from now on • denotes the usual scalar product in IR 3 and \v\ the Euclidean 
length of uglR 3 . With the abbreviation 

E:= v/l + M 2 , i. e., p° = e - fl E 

the Vlasov-Einstein system now takes the following form: 

d t f + e»- x ^-d x f-(\— +e"-VB)--a,/ = 0, (1.3) 
hi V r ) r 
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e- 2X (2r\' - 1) + 1 = 8Trr 2 p, (1.4) 
e~ 2A (2r// + l)-l = 8nr 2 p, (1.5) 



where 

p(t,r) = p{t,x):= J Ef(t,x,v)dv, (1.6) 

// x ' v \ ^ dv 
{—) f(t,x,v) — . (1.7) 

By and ' we denote the partial derivative with respect to t and r respectively. 
The phase space distribution function / is assumed to be spherically symmetric, 
i. e., f(t,Ax,Av) = f(t,x,v) for every AG SO (3) andielR, x, vGlR 3 . It should 
be noted that while the above system is well-posed is is not the complete Vlasov- 
Einstein system: Only the 00 and 11 components of the field equations are 
written. The 01 component is also used in the sequel; like the also non-trivial 
22 and 33 components it follows from the reduced system above, and it reads 

A = -47rre A+ ^, (1.8) 

where 

j(t,r)=j(t,x):= J ^f(t,x,v)dv. (1.9) 

Due to the symmetry the field equations have no radiative degrees of freedom, 
and the solution is determined by an initial condition 

f(0,X,v) = f(x,v). 

The initial data are taken to be spherically symmetric, non-negative, continu- 
ously differentiable, compactly supported, and such that 



\y\<r 



o T 

f(y,v)dvdy<-, r>0. 



The latter condition rules out trapped surfaces at t = 0. We briefly discuss the 
main analytical results on this system. In |l3[ ] it was shown that each initial 
datum as above launches a unique smooth solution for which all derivatives 
which appear in the system exist and are continuous. The solution can be 
extended in time as long as p remains bounded. For small initial data it is then 
shown that a global, geodesically complete solution results which decays to flat 
Minkowski space for t — * oc. On the other hand, it is shown in |l^] that large data 
will lead to a singularity, and in |l7| it is shown that the first such singularity 
will be at the centre. The main open problem is whether large data still lead 
to solutions which are global in Schwarzschild time; if so this would imply 
the cosmic censorship hypotheses. This problem was investigated by numerical 
simulation in pM . Given / as above the system was solved numerically with 
initial datum Af where the amplitude A > was slowly increased from small 
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to large values. There was no indication that the solutions might blow up in 
finite Schwarzschild time, and there existed a critical amplitude A* such that 
for A < A* the solution dispersed while for A > A* a black hole with mass M(A) 
seemed to form. The system exhibited exclusively the so-called type I behavior: 
\im.A—>-A,+M(A) =-M* >0. This is in sharp contrast to results for the same 
sort of numerical experiment with different matter models which exhibit type II 
behavior, M* = 0. The findings of jll| are confirmed in @] so that considerable 
evidence seems to indicate that cosmic censorship holds for the Vlasov-Einstein 
system. One should note, however, that in |27[] the formation of 

naked singularities was reported in a numerical simulation of the Vlasov-Einstein 
system with axial symmetry. On the other hand, certain numerical findings 
which were claimed there for the Vlasov-Poisson system which is the Newtonian 
limit of the Vlasov-Einstein system contradict known analytical results. 

Given the significance of these questions it is desirable to have a rigorous 
mathematical foundation for the numerical scheme which is used, that is, to 
show that the solutions of some appropriately discretized version of the system 
converge to solutions of the Vlasov-Einstein system if the parameters of the 
discretization vary appropriately, and to give corresponding error bounds. The 
aim of the present paper is to carry out this program for the scheme used in 
[ fl8| , a so-called particle-in-cell scheme. This scheme is discussed in detail in 
the next section, the basic set-up being as follows: The support of the initial 
datum is split into small cells. In each cell a point is chosen which carries a 
weight representing the integral of / over this cell. These particles are smeared 
out in space by a hat function. From this approximation for / approximations 
for the source terms in the field equations and thus approximations for the 
fields can be determined. With these we enter into the characteristic system of 
the Vlasov equation and propagate the particles by one time step. Then the 
process is repeated. This sort of numerical scheme is well known in plasma 
physics and in astrophysics where it is used to simulate the Vlasov-Maxwell or 
the Vlasov-Poisson system. For convergence results for these systems we refer 
to j|, H [| |28|, H^j] and in particular to (22) where the spherically symmetric 
Vlasov-Poisson system is considered. General background on such schemes can 
be found in @, §. 

For the present system particular difficulties arise due to the fact that the 
source terms in the field equations contribute to the fields also in a pointwise 
sense; as opposed to the Vlasov-Poisson and Vlasov-Maxwell systems they are 
not necessarily integrated in space(time). This lac k of smoothing effect of the 
field equations which is apparent for example in (l.S) causes considerable analyt- 
ical as well as numerical complications, and the analysis of a numerical scheme 
must be based on a careful analysis of the analytical properties of the solutions, 
which was initiated in pl|, cf. also p| . The paper proceeds as follows: In the 
next section we first collect some additional information on the Vlasov-Einstein 
system which is needed in the sequel. Then we state the discretized version of 
our system and the main result. At this stage we only discretize in phase space, 
thereby reducing the system to a system of ordinary differential equations for the 
particle positions and other relevant quantities. The proofs of our convergence 
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results and error estimates are carried out in Section 3 in a series of lemmas. In 
Section 4 we discuss the question of how to discretize the system also in time. 
Since our numerical findings coincide with those reported in jl7j and || we do 
not include them here. An open problem is to analyze the role which steady 
states of the system play in explaining the observed type I behavior. This will 
be the topic of a separate numerical investigation for which the present paper 
sets the theoretical stage. 

We conclude this introduction with some further references to the literature. 
In iQ and (2Cj it was shown that solutions of the Vlasov-Einstein system con- 
verge to solutions of the Vlasov-Poisson system in the Newtonian limit. The 
latter system is much better understood than the Vlasov-Einstein system, in 
particular, there exists an existence result for global smooth solutions to gen- 
eral initial data, cf. & 0, [23) . Another important feature of the Vlasov-Einstein 
system as stated above is that it possesses a large family of steady states, cf. 



[|15| |16| y_2|. All the results mentioned so far refer to the asymptotically flat 
case which is characterized by the boundary condition (1.1) and from a physics 
point of view represents an isolated system such as a galaxy or globular cluster 
in an otherwise empty universe. There exists also a growing number of results 
on the cosmological case of the Vlasov-Einstein system, and we refer to |l| for a 
discussion of and references to these results. Finally we mention that the results 
of the present paper constitute the major part of the second author's PhD thesis 



2 The semi-discretized approximation — main 
result 

o 

An initial datum / as specified above launches a continuously differentiable and 
spherically symmetric solution /. Let [0,T] denote any time interval on which 
this solution exists. For numerical investigations it is important to make use of 
the spherical symmetry in such a way as to reduce the number of independent 
variables. We introduce 

ll X ' V 2 2/ \2l |2 

r:=uc, w:= , L = x v ~{x-v) =\xxv\ 

r 

so that 

E= y/l + w 2 +L/r 2 . (2.1) 

It can be shown that / must be of the form 

f(t,x,v) = f(t,r,w,L). 

Instead of writing the Vlasov equation in these variables we write down its 
characteristic system which is equivalent but more relevant for what follows: 

r = e»~^, (2.2) 
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w = —Xw — ' 
L = 0. 



t "- x n'E + e^ x - 



r 3 E' 



(2.3) 
(2.4) 



Note that the quantity L, the modulus of angular momentum squared, is con- 
served along characteristics so that the characteristic system is essentially two 
dimensional. The field equations (1.4), (1.5), (1.8) remain unaffected by the 
above change of variables, and the source terms can be written as 



7T f°° [°° 
Pit, 7 ") = ~2 / / Ef(t,r 7 w,)dLdw, 



p(t,r) 



-oo J 

oo poo 



W 

~E 



f(t,r,w,L)dLdw, 



oo Jo 



wf(t,r,w,L)dL dw. 



(2.5) 
(2.6) 
(2.7) 



In the sequel we will use both (x,v) and (r,w,L) in our discussion, but the 
numerical scheme will be formulated in the latter coordinates. Let us denote by 
(X,V)(s,t,x,v) or (R,W,L)(s,t,r, w,L) the solution of the characteristic system 
which at time s = t takes the value (x,w)elR 6 or (r,w,L) G [0,oo[xIR x [0,oo[. 
Then (by abuse of notation) 

f(t,x,v)^f((X,V)(0 1 t 1 x,v)) = f((R,W,L)(0,t,r,w 1 L))^f(t,r,w,L). 

It is important that the field equations can be solved explicitly for A and \x and 
their relevant derivatives. Using (1.2) the field equation ( |1.4| ) can be integrated 
to yield 

-2A(t,r)_-, 2m(t,r) 



where 



m(i,r):=47r / s p(t,s)ds 



(2.8) 



(2.9) 



the right hand side of (2.£ ) is positive initially by assumption on / and remains so 
on the existence interval of the solution. It is worthwhile to note that m(t, oo) is 
a conserved quantity of the system, the ADM mass. Another conserved quantity, 
related to the conservation of the number of particles, is 



i xit - r) f(t,x lV )dvdx. 



Next we solve (|l.5[ ) for p' , 



= 2A(t,r) 



m(t,r) 



-4irrp(t,r) 



and using (1.1) this is integrated to give 



/j,'(t,s)ds. 



(2.10) 



(2.11) 



(2.12) 
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From ( 



1.4) we also have 



A'(t,r)=e 2A ^) L!^M + 47rrp(t ) r)) (2.13) 

which will become relevant shortly. Obviously, A>0, /i<0, and by adding 
( |2.11 ) and ( 2.13| ) and observing the boundary condition (1.1), X + fi<0. As 
noted above, / is constant along characteristics, but the characteristic flow is not 
volume preserving as can be seen from the factor e x in ( 2.10| ) . This complication 
is the price we pay for writing the system in non-canonical momentum variables 
which make the explicit representations of the metric quantities in terms of the 
source terms possible and greatly simplify the Vlasov equation. We need to 
know how integrals of / over pieces of phase space evolve: 

Lemma 2.1 Let Z{-,t,z) — (X,V)(-,t,x,v ) denote the solution of the character- 
istic system of ( |i.j| J with Z(t,t,z) — zElR 6 . Then 

detd z Z{s,t,z) = e x< - t ^- x ( s > x ( s > t ' z ». 

Let AclR 6 be measurable and A(t) := Z(t,0,A). Then we have for the solution 
f and any continuous function g, 

g(t,z)f(t,z)dz = f g(t,Z(tAz))f(z)e x ^ x ^°^dz 

A{t) J A 

and 

If f(t,z)dz = - f f(t,z)(\(t,r)+R(t,0,Z(0,t,z))\'(t,rj)dz. 

at JA(t) JA(t) v 7 

Proof. The first assertion follows from Liouville's Theorem, and the rest follows 
by a change of variables and straight forward computation. □ 

Since we will not be too careful in distinguishing between a symmetric subset 
AcTR 3 xR 3 and the subset A of (r,u;,L)-space which describes A, that is 

A=Ux,v)£lR?\ (\x\,x-v/\x\,\xxv\ 2 ) el}, 

we note that 

meas(A) = 47r 2 meas( J 4). 
This follows from dv = irr~ 2 dwdL and dx = Airr 2 dr . 

o 

We make the following additional assumption on the initial datum /: With 
respect to the (r,u>, L)-variables, 

o 

supp/C [r minj^max] X [w m ; n , W max ] X [L m i n , i ma x] (^--T) 

with constants < r min < r max , w min <w max , < L min < L max . Since L is con- 
served along particle trajectories, the particles remain away from zero in space 
as long as their momenta v remain bounded. This assumption, which was also 
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made in p2[ |, avoids problems arising from particles at the origin where our 
coordinates (r,w,L) are not suitable. In |H| it is shown for the Vlasov-Poisson 
system how this restriction can be avoided by switching to Cartesian coordi- 
nates in the neighborhood of the origin, but we prefer to avoid at least this 
technical complication. In passing we note that the numerical investigations in 
[ |l8| H considered data satisfying this restriction. 
We are now ready to formulate the 

o 

Semi-discretized approximation: We decompose the support of / in 



(r, w,L)-space or a supset of the form (2.14) into disjoint, connected, and mea- 
surable sets A n , n£ {1,2,..., AT}, with 

diam{A n )<e, meas(A n ) > -^e 3 , n£ {1,2,..., TV}, 

where the fineness 0<e^;l is a small parameter; the factor 1/100 can be re- 
placed by any positive constant which is then kept fixed. In each cell A n we fix 
a point (r ni w n ,L n ) £A n and define 

(R ni W ni L n )(t):=(R,W,L)(t,0,(r ni w n ,L n )), 



the corresponding solutions of the characteristic system (2.2), (2.3), (2.4). The 
cell in (x,w)-space which corresponds to A n is denoted by the same symbol. 
Moreover, we define 

A n (t):=Z(t,0,A n ), M n (t):= [ f(t,z)dz; (2.15) 

JA n {t) 



note that Lemma 2.1 tells us how M n evolves. 

Our aim is to obtain approximations to R n (t), W n (t),M n (t) as solutions of 
an appropriately set-up system of ordinary equations. Throughout the paper 
approximating quantities will be denoted with a bar. Suppose we have N points 
in (r,u>,L)-space with coordinates 



R = (i?i, . . . ,Rn) £ 1R. 



+ ' 

N 



L = (Li,...,Ln) 6lR + , 



carrying weights 

From this information we generate approximations for the source terms at any 



M=(M 1 ,...,M N )EJR, I l. 



r>0: Following (2.1) we abbreviate 



and define 



E n :=E{R n ,W n ,L n ):=^l + W 2 n + L n /R 2 n (2.16) 



i N _ 

p(r,R,W,L,M) := ^—^^£; n M nX5 (r-i?„), 

n=l 
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™™ i 



p(r,R,W,L,M) := — -J2=^M nX5 (r-R n ), 

n=l n 
1 N - 

j(r,R,W,L,M) := — — J2w n M nXs (r - R n ). 



4irr 2 5 ■ 

n=l 

Here the hat function \s which is used to smear out the particles is defined by 

x*(C)-j ) | C!>J) CeiR, 

where < <5 <C 1 is another small parameter; the necessary relations between 5 
and e will be specified shortly. Suppressing the variables R,W,L,M for the 
moment we can now define approximations for the various quantities appearing 
on the right hand side of the characteristic system etc.: 

m(r) := An / s 2 p(s)ds, 
Jo 

e -2A(r) := : _2m(r) 5 
r 

Ji(r) ■=- I ~p'(s)ds, 

J r 

\'(r) :=e 2 ^(-^-+4nrp(r)y 

A(r) := -47rre^ r ) + ^ r )j(r)- 

The approximations R n {t), W n (t), M n (t) to the true quantities R n (t), W n (t), 
M n (t) are now defined as the solutions of of the following autonomous system 
of 3iV ordinary differential equations: 

5 n = e (P-A)(fl»)^n j (2.i7) 

W n = e^- J ^^^-k(R n )W n -e^- J ^Ji\R n )E n , (2.18) 
R n E n 

i* n = - (j(Rn) + e^- J ^^\\Rn)^Mn (2.19) 

with initial conditions (R n ,W n )(0) — (r n ,w n ) and M„(0) = M„(0), ra=l,...,iV. 
Clearly, we set L n = L n - Moreover, the dependence of A, ~p, and their various 
derivatives on the coordinates and weights of all the other particles has been 
suppressed in the notation above. 

Given a solution R(t) ,W (t) , M (t) of the above discretized system it is sug- 
gestive to abuse our notation as follows: 

p(t,r):=p(r,R(t),W(t),M(t),L), 



9 



and analogously for all the other source terms, metric coefficients, and their 
various derivatives. Note that A is just a name. In general X^dtX, whereas 
~p' = d r Jl and A =d r X. Nevertheless, we will use the notation A because we will 
compare A to A during the convergence proof. 

Before we state the main result we specify certain quantities which we use to 
control the errors. By || • we denote the L°° norm with respect to r€ [0,oo[. 
For the characteristics we define 

\\R(t)-R(t)\\ :=ma X {\R n (s)-R n (s)\\ne{l,...,N}, s€[0,t}}, 
\\W(t)-W(t)\\ :=max{|W n (*)-W n (s)||ne{l,...,JV} ) se[0,t]}, 
\\M(t)-M(t)\\ :-e- 3 max{|M„( S )-M„( S )||ne{l,...,7V}, se[0,t]}; 

since the weights M n are of order 0(£r 3 ) we have inserted a suitable factor. We 
are now able to state the main result: 

Theorem 2.2 Assume in addition that f <E C 2 (H 6 ) so that the solution has the 
same regularity with respect to x and v. Then there exist constants C\ > and 
C2 > depending only on the approximated solution f restricted to the time 
interval [0,T] such that the following holds: If Q <e <5<C\ then any solution 
of the discretized system stated above exists on [0,T] and satisfies the following 
estimates for t € [0,T]: 

\\R(t)-R(t)\\ + \\Wit)-W{t)\\ + \\M(t)-M(t)\\ <<?(*» + £), 
|| A(t) - A(t)|| 00 + \W)-m\\oo + \\m(t) - ro(t)||oo < C (s 2 + 1) , 

\\p(t)-m\\o +\\p(t)-p(t)\\ 00 +h(t)-m\\oo<c{6+^). 

Note that the error estimates for the source terms are one order worse than the 
other estimates; this is no coincidence. 

Remark. Strictly speaking there exists in the literature no proof for the 

o o fi 

assertion that the regularity assumption feC (H ) will propagate to the so- 
lution. However, d Xi f will satisfy a system of differential equations which is 
obtained from the Vlasov-Einstein system by differentiating the various equa- 
tions accordingly, and this system will be linear in d Xi f so that this derivative 
will be C 1 on the existence interval of /. Turning this into a proof would be 
a lengthy exercise. A version of our theorem without the additional regularity 
assumption is also possible. However, we then no longer obtain error bounds for 
the source terms, and instead of the regularity assumption an a-priori bound of 
the approximating source terms, uniform in e and 5 has to be made which can 
be monitored during a run-time situation. Since the statement of this second, 
weaker form of our result requires more technical preparation we postpone it to 
the next section. 
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3 The semi-discretized approximation — proofs 



First we need to introduce some further information on the solutions of the 
Vlasov-Einstein system. From jl3| we know that any initial datum JeC 1 (]R 6 ) 
as specified above launches a unique solution /. Let [0,T] be any time interval 
on which this solution exists. Then /GC 1 ([0,T] xIR 6 ) with respect to (x,v) 
and 

ftP,jeC 1 ([0,T] x [0,ooD, m,X^eC 2 ([0,T] x [0,oo[) 

as functions of t and r. This solution together with the time interval are now 
kept fixed, and we will need the following bounds: 

Lemma 3.1 There is a constant D>1 such that for all t 6 [0,T] and r>0, 

p(t,r), p(t,r), \j(t,r)\, e 2A ^\ |A(t,r)|, |A'(*,r)|, 

\»(t,r)\, \n'(t,r)\, \X'(t,r)\, \X"(t,r)\, |^"(t,r)| < D, 

o 

and for all t £ [0,T] and (r,w,L) € supp/, 

R{t,0,r,w,L), — — -, \R(t,0,r,w,L)l \W(t,0,r,w,L)\<D. 

R{t,0,r,w,L) 

Moreover, D can be chosen such that 

e- 3 M n (t)<D, n=l,...,N, 



for any discretization o/supp/ as specified above. 



Proof. The right hand side of (2.2) is bounded by 1, and hence the source terms 
vanish outside the light cone, that is for all (t,r) with r>r max + t. Together 
with the regularity of the solution this implies the first set of estimates. Since 

o 

V(t,0,z) is uniformly bounded for t G [0,T] and z esupp/ and L is conserved we 



get a lower bound on R. Finally, by Lemma 2.1 



M n {t)= f(t,z)dz< f(z)e x(0 ^dz<\\fe x ^\\ oc me aS (A n ), 

JA„(t) JA n 

and the proof is complete. □ 
All the functions on the right hand side of (2.17), (2.18), (2.19) are Lipschitz- 
continuous in the unknowns. Let [0,2^ denote the maximal interval on which 
the approximate solution exists and satisfies the bounds 

e- 3 M„(t),i?„(t),=--,|TF„W|,||e 2 ^|| 00 <2 J D, n=l,...,N. (3.1) 
R n (t) 



Here D is the constant used in Lemma 3.1. It must be emphasized that T Ej s 

o 

depends on the particular discretization of the support of / and not just on e 
and 5. 

We can now formulate the weaker version of our main result which does not 
require the additional regularity of the approximated solution: 
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Theorem 3.2 Let a family of approximating solutions parameterized by e and 
5 (or a sequence of such) be given and assume that for some constant C* > 0, 

- P {t)<C\ 0<£<<5<^, ie[0,T e ,,]. 

Then there exist constants Ci>0 and C2 > depending on the restriction of 
the approximated solution f to the time interval [0,T] and on C* such that 
the following holds: If < e < S < C\ then any solution of the discretized system 
stated above exists on [0,T] and satisfies the following estimates for t G [0,T]: 

\\R(t)-R(t)\\ + \\W(t)-W(t)\\ + \\M(t)-M(t)\\ < c(s+ £ - 

||A(*) - A(t) Hoc + H^(t) -p(t) || oo + l|m(*) -5IS(t) Hoc <C (<* + !)■ 

Note that in a run-time situation the assumption on the boundedness of ~p can 
be monitored by looking at the numerical data. 

Throughout the rest of this section we fix discretization parameters 

0<e<6<—. 
~ ~ AD 

The lower bound on vneas{A n ), n=l,...,N, which is to hold for any discretiza- 
tion implies that 

Ne 3 <C; (3.2) 

as in all the estimates which follow C denotes a positive constant which may 
depend on D, C* , or the restriction of / on [0,T]. We assume that the assump- 
tions of Theorem |2.2| or Theorem 3.2 hold. In the situation of Theorem 2.2 we 



redefine T e j such that [0,T £i( $] is the maximal interval on which the estimate 

||p(*)lloo<2£> (3.3) 
holds in addition to those stated in (13 . If) . In the arguments which follow there is 



almost no need to distinguish between Theorem 2.2 and Theorem 3.2. The o nly 
difference is the bound on p: p(t, r) < 2D for 1 6 [0,71 5] in case of Theorem |2.2| 
and ~p(t,r) < C* for t G [0,T ej ,s] in case of Theorem 3.2 . 

We will use the following abbreviation for the index set of the decomposition 

o 

of supp/: 

We start our chain of auxiliary results by collecting some information on our 
hat function \s'- 

Lemma 3.3 Let i5>0 and define 

rC 
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Then 

|x(C)-x(OI<JlC-£l, 0<x(C)<i, C, eeM, 

and for a,b,C G El wif/i |a — C| > \a — b\ + 5, 

x(C-a) = x(C-6). 

The proof is straight forward. Next we collect some additional information 
on the characteristic flow of the Vlasov-Einstein system. As before, we often 
abbreviate z = (x,v) € IR 3 x IR 3 and Z — (X, V) accordingly. 

o 

Lemma 3.4 For all characteristics starting in the support of f and all t £ [0,T], 
\Z(t,0,z 1 )-Z(t,0,z 2 )\<C\z 1 -z 2 \, 

O 

in particular, for A C supp/ C IR and A(t):= Z(t,0,A), 

diam(A(t)) < Cdiam(A). 

Also 

\{R,W){t,Q,rx,w x ,L x )-{R,W){t,Q,r 2 ,W2,L 2 )\ 

<C(\r 1 -r 2 \ + \w 1 -w 2 \ + \L 1 - L 2 \). 

Proof. The first estimate is included in O, Prop. 2.2], with the bounds from 



Lemma 3.1 this implies the third estimate; note that 1/R{t) is bounded. □ 
We now start establishing bounds on the approximating solutions; recall that 
constants C may depends on the true solution /, but never on the discretization 
parameters e and 5. 

Lemma 3.5 For all 1 6 [0,T ejl j], r > and neJ\f, 



p(t,r),p(t,r), \j{t,r)\,m(t,r),\X(t,r)\, \p{t,r)\, \X(t,r)\, |A (t,r)\, \jjf(t,r)\<C 



M n (t) 



<Ce i 



Moreover, for r [1/(4D),3D], 

p(t,r)=p(t,r)=j(t,r) = 0. 

Proof. Let te[0,T £iS ]. By definition of T s j, r<l/{AD) implies r-R n (t)< 
l/(4L>)-l/(2L>)<-'<5, and r>3L» implies r - R n (t) > 3D - 2D > S. Thus 

p(t,r) = 0, r^[l/(4D),3U]. (3.4) 

A bound on ~p of the desired type holds either by assumption or by definition of 
T e ,s- It is obvious that p, \j\ < ~p so that the assertio ns fo r p and j follow. Using 



the bounds on the approximating source terms and ( ^4 ) the remaining bounds 
follow from the definitions of the respective quantities and of T e ^- □ 
Next we present an important tool for analyzing the approximating source 
terms: 
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Lemma 3.6 LetO<r <ri. For t e [0 ,T E t s) define 

I(t) :={neJN\Rn(t) G [r ,n]}, J(t) := {n6H|^,(t)e [r ,n]}. 

TTien 

|/(*)| < cn^f±f , 

|7(t)| < c r 1 -r + p(t)-^)||+ £ 

where \I\ denotes the number of elements of a set I cIN. Moreover, 

J2 M n (t)<C(n-r + 5). 

nel(t) 

Proof. The proof of the first assertion follo ws the one of j(| Lemma 2]. Fix 



n£l(t) and (r,w,L) &A n . Then Lemma 3.4 implies that 
R(t, 0,r,w,L)\ £ [r -Ce,ri + Ce], 

and thus 

(J A n {t) c [r -Ce,n+Ce} x [-D,D] x [0,D] 



Using Lemma 2.1, the bound on A from Lemma 3.1, and the lower bound on 



meas(A„) we find 

{r l ~r + 2Ce)2D 2 > ^ meas(A n (») > - ^ meas( J 4„) > C\I{t)\e 3 , 

nel(t) n€l(t) 

which implies the first assertion. The second one follows from the first, since 

T(i)c{nGA/'|i2 n (i)G[ro-||i2(i)-S(t)||,ri + ||fl(t)-S(t)||]}. 

To prove the third inequality define for r > 8/2, 

7(t,r) := {neJV| R n (t) G [r-£/V + <V2]}. 

Then nel(t,r) implies |E„(t) - r| < 5/2 and therefore xs(R~n(t) -r) > 1/2. 
Hence 

5^ M n (t)<2 M n (t)xs(Rn(t)-r) 

< 8Trr 2 5p(t,r) <8vrr 2 5m&x{C* ,2D} = Cr 2 6, r>5/2. 

Now define k :— [(r x —r )/S~\ and :=To + i£, i = 0, ...,k. Here [CI denotes the 
smallest integer larger than or equal to (. Then sq = ro and Sk > v\, and therefore 

fe-i 



I(t)c\Jl{t, Sl + 5/2). 



i=0 
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Using this inclusion and the previous estimate for every i together with the 
fact that by definition of T Ej s only r\ < 3D needs to be considered we get the 
estimate 

fc-i fc-i 
^M„W<]T ]T M n (t)<J2C( Sl + 5/2) 2 5 

neT(t) i=Q neI(t,Si+S/2) i=0 

< C(n + 8/2) 2 kS+ < C |"(n -r )/<5]5< C(n - r + <f) , 

and the proof is complete. □ 

In order to analyze the errors appearing at the source terms and the metric 
coefficients due to the discretization we define quantities that are intermediate 
between the true quantities and their approximations in the sense that the 
formulas for the latter source terms are evaluated at the true characteristics 
with the true weight functions. These intermediate quantities are denoted by 
double-bars: 

1 N 



4:Ttr 2 5 ■ 

n—l 



where analogously to fl2.16| ), 



E n (t) := y/l + W*(t) + L n /H*(t), 

with corresponding definitions for p and j. Moreover, M n is to be the solution 
of 

H„ (t) = - (X(t, R n (i)) + R n (t) A' (*, (*))) H„ (t) , W n (0) = M n (0). 

Note that we have, similar to p and M, \\p(t) +e" 3 M(t) < C for all te[0,T] 
and^(t,r)=0 for all r£[l/(4D),3D], te[0,T]. 

Lemma 3.7 For all 1 6 [0,T e . s ], 

|IM*)-7*(*)l|oo + IIAW-ACOHoo + llmW-mWHoo 
<c(||/j(t)-^(t)||oo + ||p(t)-f(t)||oo 

+ \\R(t)-R(t)\\ + \\W(t) -W(t)\\ + \\M(t)-M(t)\\ 

Proof. Let t G [0,T £j i] and r > 0. We start with analyzing the error m — m. Due 
to the definition of m, /?, and p, 

m(t,r)-m(t,r) = 4:TT / (p(*,s) - p{t,s))s 2 ds = 4tt / (p(t,s) -p z (t,s))s 2 ds 
Jo Jo 

i r N 

+ S J I]( S ™W^r l WX5(s-i?nW)-^nW^r l (t)X5(s-i?„(t)))d S 

=:Fi+F 2 . 
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Since for r > 3D we have p(t,r) = p(t,r) = 0, 

|*i|<c||p(t) -?(*)!!«, 

Using the definition of x, cf. Lemma |3.3|, we have 



IF, I < 



N 



53 (M„ (t) - M n {t))E n (t) X (r - Rn (t) ) 
i 

N 

Y^n{t) (E n (t)-E n (t)) X (r-Rn(t)) 

71=1 

AT 

53 M n (t)E n (i) ( X (r - - X (r - 1„ (t))) 



! F21 + -F22 ' 



Since E n and x are bounded and e 3 N<C by (3.2) 



F 21 <Cj2\M n {t)-M n (t)\<C\\M(t)-M(t)\\. 

n=l 

Since W n (t), W n (t), 1/R n (t), 1/R n (t), and L n are bounded, 

I (t) - E n (t) I < C (I iin (t) - R n (t) I + 1 W n (t) - W„ (i) I) . 

and thus 

JV 

f 22 < c53M n (t)(iiij(t)-s(t)ii+nw(t)-w(t)ii) 

n=l 

<C(\\R(t)-R(t)\\ + \\W(t)-W(t)\\). 

If we define 

7i(t):= {ne^|| J R n (t)-r|<]]i?(t)-l(t)|| 00 + <5} 
then by Lemma |3.3| , 

X (r - i?„ (t ) ) = x (r - R n (t ) ) , n & h (t ) , 



and by Lemma 3.£ 



|A(*)|<C 



which implies, again with Lemma |3.3|, that 



F 23 <C ^ M n (t)min\ lRn{t) Rn{t)l ,l\ <C\\R(t) 



-R(t)\\. 



neii(*) 
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This completes the desired estimate for m — m. The estimate for A — A follows 
directly from the expressions for A and A in terms of m and m, the fact that 
A(t,r) = X(t,r) = for r < 1/(4D), and the estimate which we just completed. 

It remains to consider /i — 71. For r < 1/(4D) we have n'{t,r) = Jl'(t,r) — and 
therefore \fi(t,r)-Ji(t,r)\ = \iJ,(t,l/(AD))--p(t,l/(4D))\. Let r > 1/(41?). Then 
the definitions of /z, 71, /it', and 7i' yield 



|/z(*,r)-/i(i,r)| 



fj, (t,s) — /i (t,s)ds 



< 



a 2X(t.s) 



m(t,s) m(t,s) 



,2\{t r 



_ _2A(t,s) 



fi.S 



m(t,s) , , 

3 ; +47TS P ft,s) )ds 



47rse 2A(t ' s) (p(t,s)-p(t,s))ds 

r //iW-mWHoo + CllAftJ-AWHoo 

47rse 2A(M) (p(t, s) -p(i, s)) ds 



+ 



We split the difference p — p into p — p and p — p. The integral resulting from the 
second term can obviously be estimated as desired. As to the first term, 



Anse 2X ^ s) (p(t,s)~p(t,s))ds 

r 3D 



< J2\Mn(t)-M n (t)\ \ 

6 n =l Jl/(4D) 



Xs(s-R n (t))ds 



c 



En(t) E n {t) 



3D 



1/(413) 



X5(s-i?„(t))c?S 



n— 1 
:F 3 + F 4 + F 5 



oo e 2A(i,s) 



{Xs(s - R n (*)) - * 4 (a - #n (*))) 



where we used the bounds on the supports of p and p and recalled that 
f>l/(iD). The first two terms are comparati vely easy to deal with: De- 
fine I 2 (t,s):={neAf\ \s-R n (t)\<6}. By Lemma |J, \I 2 (t,s)\ < C(S + e)e- 3 < 
C5e~ 3 for all s>r. Thus 



F 3 <je 3 \\M(t)-M(t)\\ 



■3D 



l/(4£>) 



|I a (t,«)|ds<C7||M(t)-Af(t)||, 



F 4 < 



1 (|| - i2(t)|| + 1| - W(*)||) 



3D 



l/(4£>) 



\h(t,s)\ds 



<C(\\R(t)-R(t)\\ + \\W(t)-W(t)\\). 
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To deal with F$ we fix some n£j\f and assume without loss of generality that 
Rn(t) < R n (t). A change of variables yields 



oo 2A(M) 



(XS (S - Rn (*)) - XS (S - Rn (*))) ds 



< 



r--Rn(t) Z + 
r-R„(t) e 2A(t,z+ii„(i)) 

r-R n {t) Z + R n (t) 



Xs(z)dz- 
■Xs(z)dz 



a 2\(t,z+R n (t)) 



-il„(i) Z + -Rn(*) 



-Xs{z)dz 



<C 



r-R n {t) 
r-~R n (t) 

r-R n (t) 



a 2\(t,z+R n (t)) 2\(t,z+R. n (t)) 



Z + R n (t) 



Xs(z)dz + C 



Z + Rn{t) 



Xs(z)dz 



r-R„(t) 



Xs(z)dz\R n (t)-R n (t)\ 



Define 



< Cmin { \R n (t) - R n (t) \,5}+ C5\\R(t) - R(t) \\. 

I 3 {t) := {neM\[r-R n (t),r-R n (t)}n[-S,+5}^9} 
C {neAf\ \R n (t)-r\<S+\\R{t)-R(t)\\}. 



Then Lemma |J yields \h(t)\ < Ce~ 3 {5+ \\R(t) -R(t)\\), and putting every- 
thing together, we obtain 

(7 F 3 _ f F 3 N _ 

1^1 <— E ™n{\R n (t)~R n (t)\,S} + — Y / S\\R(t)-R(t)\\ 

nel 3 (t) n=l 

<C\\R(t)-R(t)\\ 

by (|3.2|), and the proof is complete. □ 

Next we estimate the differences of the intermediate, double-barred source 
terms and their true counterparts: 

Lemma 3.8 For all t G [0 ,T e i<5 ] , 

||p(0-^(f)lloo <C(S + e/S), 

||p(*)-?(*)||oo < C(S 2 +s/5), provided p"eC([0,T]xE + ), 
and the analogous assertions hold for p and j. 

Proof. We restrict ourselves to considering p; p and j can be dealt with in the 
same fashion. Let t€ [0,T ei g] and r>0. Then 



4irr 2 5 p(t,r) 



r+8 



Airs 2 p(t, s)xs(r - s) ds 



-5 
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r+S 

+ j 4:Tr(r 2 p(t 7 r)-s 2 p(t,s))xs(r-s)ds 



Ir-S 

N 



4nr 2 Sp(t,r) = Y^E^xs^ - R n {t)) j f(t,z)dz 
„=i J A n {t) 

N 

= V / Ef(t,z) X 6(r-\x\)dz 

n=1 JA n {t) 
N 

+ V / (E n (t)-E)f(t,z)xs(r-R n (t))dz 
n=1 JA n (t) 

N 

V / B/(t, 2 ;)(x S (r-i2 n (t))-X«(r-|a;|))<fa 
rr; Jam 



= : ^3 + ^4 + ^5. 

Clearly, Fi=F 3 . Since p(t,r) = £(t,r) =0 for all r [1/(4D),3D], 
|K*,r)-p(t,r)|<^(|J!,| + |fi| + |F 6 |). 

Using the fact that we assume a bound on the derivative of s 2 p(t, s) with respect 
to s we find 

rr+S 



\F 2 \<C f \r- S \xs(r-s)ds<CS 2 

Jr-S 



To estimate F4 let z— (x,v) £A n (i). Let = (X n (t),V n (t)) be any Cartesian 
characteristic starting at a point with (r,w,L) coordinates (R n (0),W n (0),L n ). 
Then Z n (t)eA n (t) and \z - Z n (t)\ < diam(A„(i)) < Cs by Lemma [3[ The 
function E used to compute p and p is Lipschitz on the relevant domain; 
the same is true for p and j. If we define Ii(t) :—{n^M\ \R n (t) — r\ <5} then 
Lemma 3.6 yields <C5e~ 3 and therefore, 



raG/i(t) 

< 



n£l l (t) 

For F5 we define 



F 4 \ < f \E n (t)-E\f(t,z)dz 

V Cdiam(A„(i)) / f{t,z)dz<C5e. 



I 2 (t):={n^M\ \R n (t)-r \ < 8+ maxdiam(A m (i)) 

^ m6-<V 

Then xs( r ~ Rn(t)) — Xs( r ~ \ x \) f° r au n^I^if) and zeA n ((), and by 
Lemma 3.6, \l2{t)\<C£~ 3 (S + m&x me _\fdi&m(A m (t))+e)<Ce~ 3 5. Hence 

1^1 ^lE f Ef(t,z)\R n (t)-\x\\dz 
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< 



c 



dmm(A n {t))M n {t)<Ce, 



n€l 2 (t) 



and the first assertion is proven. 

To prove the second one, we have to sharpen the estimate for F2 in the 
case that p" exists and is continuous and thus bounded on the relevant domain. 
First we rewrite this term; since throughout the argument t remains fixed we 
suppress this variable: 

rr+8 

F 2 = 4tt / (r 2 - s 2 )p(r)xs(r- s)ds 

Jr-S 

+ 4tt^ (p(r)-)^(p{r + 5)+p{r-8) S j S j s 2 X 8{r-s)ds 
\{p{r + 5)+p{r-5)^j -p{s)j s 2 X s(r-s) 



-4-7T 



r+8 



r-8 



ds 



='■ Fe + Ff + Fg. 
By the mean value theorem 

\F.\<C 



r+8 



r — S 



(r 2 ~s 2 )ds 



■-C5 3 . 



To analyze Fy we use the extra regularity of p. By Taylor expansion, 



\F 7 \<cj^ P (r)-\[p(r + 5)+p(r-5)) 



s 2 ds<C5 3 . 



Finally, by the mean value theorem and Taylor expansion, 

f-r+8 



\F 8 \ < C 



C 



■+s 



< C 



J s (^(p{r + 6) + p(r-Sfj-p(s)^ds 
Sp(r + S) + Sp(r-S)- 

rr+8 

Sp(r + S) + 5p(r-5)- / (p(r + S) + p'(r + 5)(s-r-6))ds 

J r 



p(s)ds— I p(s)ds 

'r-8 



r — 8 



(p(r — 6) + p (r — S)(s — r + 5)) ds 



rr+8 rr 

+ C / (s-r-5) 2 ds + C / (s-r + S) 

J r J r — 8 



'ds 



= C 
< C8 3 



l -5 2 p'(r + 5)- l -5 2 p<(r-5) 



-CS 3 
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and the proof is complete. □ 

The next lemma will complement the previous one in that we now estimate 
the differences between the approximations of the source terms and the cor- 
responding intermediate double-barred quantities. It is at this point that we 
loose one order of 5 in the error estimates for the sources compared to the other 
quantities in Theorem 2.2: 

Lemma 3.9 For all t G [0 ,T Si g\ , 

WW) -m\\«, < c OMh-^M + \\w(t)-W(t)\\ + \\M(t)-M(t)\\ 

in particular, ||p(0) — /o(0) ||oo = 0. The same estimates hold for p and j. 

Proof. Since ~p(t,r) = p(t,r) =0 for r ^[1/ (AD), 3D] we only need to consider 
re [I /(AD), 3D]. By definition 

1 N _ _ - 

p(t,r)-p(t,r) = —^J2 E n(t)M n (t) (xs(r-R n (t))-Xs(r-R n (t))) 

71=1 

1 - - — 

+ Y ( En w m " (*) - En (*) M » (*)) xs ( r - Rn (*) ) ■ 

n=l 

Let h(t):={neJV\\R n (t)-r\<5} and I 2 (t):= {neM\\R n (t)-r\<S}. Then 
Lemma |3^ implies that \h (t)\ < Cs~ 3 S and J2nei 2 (t) M n (t) < C8. Moreover, n $ 
I\ (t) U (t ) implies xs (r — Rn (t) ) = = \s ( r — Rn (t ) ) ■ Using similar estimates as 
in the lemmas above we find 

\-p(t,r)-=p(t,r)\<^ Y, Mn(t)\Rn(t)-R n (t)\/8 
neli(t)ul 2 (t) 

(J £ 3 __ 

+ — Y (\Rn(t)-Rn(t)\ + \W n (t)-W n (t)\+^ 3 \Mn(t)-M n (t)\) 
neh(t) 

< c \\R(t)-m\\ ( E £3+ E 

\neli(t) n£l 2 (t) 

r< F 3 

+ — \h(t)\ (\\R(t)-R(t)\\ + \\W(t)-W(t)\\ + \\M(t)-M(t)\\) , 

and the proof is complete. □ 
We now turn to the analysis of the evolution equations. 

Lemma 3.10 For all te [0,T e ,s\, 

\\M(t)-W(t)\\ :=£- 3 max sup \M n (s) -H„(s)| <Ce. 
neA ^ se[o,t] 
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Proof. Let t<E[0,T e ^] and ndAT. The definition of M n and Lemma 2.1 with 
A = A n yields 

M„(*)-If n (i) = f ^X(s,R n (s))+R n {s)X'(s,R n (s)fjW n (s)ds 



f(s,z)(X(s,\x\) + R(s,0,Z(0,s,z))X'(s,\x\)) dzds 



ds 



J A n (s) 

(x(s,R n (s)) + R n (s)X' (W n (s)-M n (s)^ 

f(s,z) \x(s,R n (s)) - X(s, \x\) + R n (s)X'(s,R n (s)) 

'0 JA n (s) 1 

- R(s,0,Z(0,s,z))X'(s,\x\) dzds. 

For z = (x,v) ei„(s), 

|A(s,i?„(s)) - A(s, |x|)| < C*|i?„(s) - |x| | < Cdiam(yl„(s)) < C*e, 

and 

1^(5)^(5,^(5)) - R(s,0,Z(0,s,z))X'(s,\x\)\ 

< R(s,0,z n )-R(s,0,Z(0,s,z))\\X'(s,R n (s))\ 

+ \R(s,0,Z(0,s,z))\\X'(s,R n (s))-X'(s,\x\)\ 

< C(\z n -Z{0,s,z)\ + \R n (s)-\x\\) 

< C(diam(A„(0)) + diam(A„(s))) <Ce. 

Hence, 



M n (t)-M n (t) 



< C 



< C 



M n (s)-M n (s) ds + Ce / M n (s)ds 



M n (s)-M n (s) 



ds + Ce 4 , 



and a Gronwall argument completes the proof. 
Lemma 3.11 For all t e[0,T e j] andneAT, 

\R n (t)~Tl n {t)\ < c(\R n (t)-Rn{t)\ + \W n {t)-W n {t)\ 
+ [[/i(t)-7Z(f)||oo + ||A(t)-A(t)|| <> 



□ 



Proof. By the equations for R n (i) and R n (t) and straight forward estimates, 
\Rn(t)-%{t)\ < 



e ( M -A)(i,fl n (t))_ e (/*-A)(t,R n (t)) 



e (M-A)(t,fl„(t)) _ e (7I-A)(t,fl„(t)) 



W n (t) W„(t) 
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and the assertion follows. □ 
The use of the following technical lemma will become obvious later: 

Lemma 3.12 Let g n<m e C l ([Q : T eiS \), m,neJ\f. Then for all t e[0,T £t s] and 
me AT, 

l N ft r _ „ v 



- X s(R m (s) -Rn(s)) [ R m {s)-R n (s) 
< Ce' 3 max (\\g k>l \\ 00 + \\g k j\\oo) \\R(t) -R(t)\\. 



ds 



fc.ieJV 

(Note that C does not depend on the functions g m ,n-) 

Proof. Let to 6 TV be fixed, denote the left hand side of the inequality by S, 
and define 

d n (s):=R m (s)-R n (s), d n (s):=R m (s)-R n (s), a€[0,T e>s ]. 
Integration by parts and d n (0) = d n (0) implies that 

S =Y1 9n,m(t)(x(d n {t))-x(d n (t))) - / ffn,m(s)(x(rfn(s))-X( d n(s))) rfs 
n=l L JO 

For se[0,T £)5 ] define 

7(s):={neAr|| J R n (s)-l m (s)|<5+||fl(s)-S(s)||}. 



By Lemma 3.6 



\I(s) | < Ce- 3 (5+ \\R(s) - R(s)\\+e) < Ce~ 3 (S + \\R(t) - R(t)\\) , s e [0,T e , 5 ] . 

If R n (s)<R m (s) — 5— \\R(s) — R(s) || then d„(s)><5 and rf„(s)><5, and hence 
x(d„Js)) = l = x(^(>)). If i?„(s)>i? m (s)+_5+||i 1 , ( S )-i?( S )|| then d„(s)<-c5 
and d„(s) < —S, and hence x(^n( s )) = = x(^n( s ))- Hence, using Lemma [D| 



\S\ < \9n,m(t)\\x(d n (t))- X (d n (t))\ 

nel(t) 

- max|3 fciZ (s)| V |x(d n (s)) -x(d»(s))| ds 



< max ||<7 M |/(t) | mm{ ||JZ(i) - S(t) || A 1 } 

+ max ||<7 M || 00 / \I(s)\mm{\\R(s)-R(s)\\/S,l}ds 
k,ieN J 

<Cm|x(|| Sw || 00 + ||sfc,i||oo) (6+\\R(t)-R(t)\\)e- 3 
mm{\\R(t)-R(t)\\/5,l}, 
and the assertion follows. □ 
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Lemma 3.13 Let 



Then 



eO) := ||p(*)-?(*)J«, + |b(a)-p(a)lloo + llj(a)-3(*)lloo 
+ ||W(s) -F(s)|| + \\M(s)-M(s)\\. 

\M{t)-M{t)\\<C f e(s)ds + C\\R(t)-R(t)\\, te[0,T e , s ] 



Proof. Let t€ [0,T e< s] and me Af. By the definitions of M m and M m , 

W m (t)-M m (t) 

A(s,S m (s))+7? m (s)A'(s,i? m (s)) M m (s)ds 



X(s,R m (s))+R m (s)X'(s,R m (s))^M m (s)ds 

A(s,l m (s))+i m (s)A'(s,S m (s)) 

-A(s,i? m (s))-i? m (s)A'(s,i? m (s)) M m (s)ds 
X(s,R m (s))-X{s,R m (s)) + (^(s)-R m (s)^X\s,R m (s^ 
+E TO (s)(A / (s,S TO (s))-A'(s !J R m (s)))lM m (s)ds 



+ 



A( S ,i? m (s)) + i? m (s)A'( S ,i? m ( S ))] (M m (s)-M m (s) 



ds 



By Lemma 3.1 



|^ 3 |<C / M m (s)-M m (a) ds. 
Jo 



By Lemma 3.1 and Lemma 3.11 

\F 2 \ < Ce 3 f (\R m (s)-R m (s)\ + \R m (s)-R\ n (s)\)ds 
Jo 

< Ce 3 ^ (\M*)-n*)\\co + \\K*)-*(«)\\°o 

+\\R(s) -R(s)\\ + \\W(s)-W(s)\\) ds. 

Now we use the formulas for the derivatives of A, A respectively where we smug 
gle in the intermediate double-barred source terms to estimate F^. 



4tt / R m (s)M m (s) 
o 
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< Ce~ 



e^ + ^ s ' Rm ^j{s,R m {s))-e^+^ s ' Rm ^j(s,R m (s)) 



2\(s,R m (s)) rn(s,R m (s)) 2 \(s,R m (s)) m { s ,Rm(s)) 



ds 



ds 



+ 



4tt^ i? m (s)M m (s)[( 



,(A+/x)( s ,fl m ( s )) _ (A+77)( s ,fl m (s)) 



)j(s,R m (s)) 



-Rm(s) (t 



2\(s,R m (s)) _ 2\(s,R m (s)) 



^jp(s,R m (s)) 



ds 



+ 



4tt/ R m (s)M m (s) 
Jo 



(A +7 i)( s ,fl m ( s )) ( j(Sj i? m(s)) _^ (s ^ ro(s) )) 

4tt / i? m ( S )M m ( S )[ e ^< s ^( s »0( S ,fi m ( S ))- J ( S ,i? ? 4 S ))) 
-7? m ( S ) e 2 ^ s ^( s ))(^( S ,i? m ( S ))-p(,s,i? m ( S ))) 



The troublesome term is the last one containing the differences j — j and p — p. 
If we denote it by F 4 , then 

l^i I < Ce 3 [ f||m(«)-?n(a)|| 00 + ||/i(«)-7l(«)||oo + ||A(s)-A(«)|| 00 

+ ||p(*)-p(*)||oo + lU(s) -!(*)||oo] ^+ \Fi\. 

To continue we introduce some abbreviations: 

p 2X(s,R m (s)) _ _ 

h m (s):= - - M m (s), ffn ( S ):= e ^- A >( s ^ s », s n ( s ):=e<"- A ><*>*««>. 

The definitions of and 7^ are analogous to those of g n and g n with J?„ instead 
of Rn, and in particular 

W n (s) ^ W n ( S ) _ ^ 

:9n(s) = Rn(s), ■=- TT g-{s) = R n {s). 



E n {s) ' 



E n {s) ■ 



25 



Moreover, H^ n (s) ~ X s(R m {s) - R n {s)). The definitions of H^ n , H^ -, and 
H^ m — should be obvious. Inserting these abbreviations and the definition of the 
source terms into F4 yields 



1*41 = 



] l E M*) (w n M n H^ n - W n M n H^) (s 

-± m (s) {E n M n H^ n -E n M n HL^ ( 
I ft N _ _ 

I ft N 

T / J2 h "^9w(s)Hi, n (s)(W n M n -W n M n )( S )d S 

° J n=l 

1 ft N , ~E \ 

- / ^h m (s)M n (s)H^ n (s)(g n W n -g n W n ^\(s)ds 

JO n=1 \ n J 

t N _ 

h m {s)R m {s)H^ n {s) (E n M n - E n M n ) (s) ds 

n=l 

1 ft *r __ _ 

1 r* N - — 

7 / Yh m ( S )E n {s)M n {s) 

[HL n ( s ) (R n - t m ) ( s ) - HL-(s) (t n ~t m )( s ))d s 

—■ F5+F6+F7 + F8 + F9 + Fio . 

Let h(s) := {neMUR n {s)-R m (s)\ < 6}, I 2 (s) := {neAf\ \R n (s) ~R m (s)\ < 6}. 
Then by Lemma |j, \h(s)\ < Ce~ 3 S and E nG / 2 ( s ) M n (s) < C5. Hence, 



\F 5 \ <C £ - f Y \R n ( s )-R n {s)\d s <Ce 3 f \\R(s) -R(s)\\ ds, 
" Jo TT 1 , , Jo 

\F 6 \<CjJ^ £ (\W n (s)-W n (s)\+e- 3 \M n (s)-M n ( S )\y s 



neli(s) 



< Ce A J [\\W(s)-W(s)\\ + \\M(s)-M(s)\\jds, 

^- C jf E (|IMs)-m(s)IU+||a( s )-a( s )|| 



n6Ji(s) 



+ \R n (s)-R n (s)\ + \W n (s)-W n (s)\)ds 
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<Ce 3 J (||m(s)-m(s)[[oo + ||A( S )-A( S ) 



\\R{s)-R{s)\\ + \\W{s)-W{s)\\)ds, 



\Fi\ <C- V (\Rn(s)-Rn(s)\ + \Wn(s)-W n (s)\ 

+e- 3 \M n ( S )-M n ( S )\^ds 
< Ce 3 j' (\\R(s)-R(s)\\ + \\W(s)-W(s)\\ + \\M(s) -M(s)\\) ds, 

£ 3 ft _____ 

\ F v\< C J Yl M n (s)\R m {s)-R n (s)\ 

mm{\R n (s)-R n (s)\/S,l}ds 
^ C jf £ M n (s)(S+\\R(s)-R(s)\\) 

mm{\\R(s)-R(s)\\/6,l}ds 
<CjJ\\R{a)-R\e)\\ _>__ M n (s)ds 



neI 1 (s)Ul 2 (s) 



< Ce 3 / \\R(s)-R(s)\\ds. 



To deal with F w let g n , m (s) :=h m (s)E n (s)M n (s). By Lemma ___| \M n \ <Ce 3 , 
hence for all n^mEAf, 

||^M,m || oo 

<Ce 6 , \\g n , m \\oo<Ce 6 . 



Thus Lemma ft.l2| yields 

|*io| <Ce- 3 mM(|| flfcl ,|| 00 + ||fl fclI || 00 )||Ji(t)-S(t)|| 
< Ce 3 \\R{t)~R(t)\\. 



Putting all our estimates together and observing Lemma 3/7, we get 



\M m (t)-M m (t)U c ^|M m ( S )-M m ( S )| rfg + c r e{s)ds + cm) _ mi 



-/o £ 3 Jo 



Taking the maximum over m G A/" and using a Gronwall argument completes the 
proof. □ 



Lemma 3.14 With e(s) defined as in Lemma 3.1\ 



\\W(t)-W(t)\\<C e(s)ds + C\\R(t)-R(t)\\, te[0,T e ,«]. 
Jo 
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Proof. The proof is similar to the one of Lemma |3.13| . Let t<E[0,T ei s] and 
me J\f. From the differential equations for W m and W m , 



W m (t)-W m (t) 



(p-\)(s,R m (s)) T 

RUs)E m (s) 



-\{s,R m (s))W m (s) 



_{JE-A)(»,fl m (a)) T 

-eMW)V( A (,))4(,)- e _ 3 _ ^ 

E ro (s)F m (s) 

+I( S ,i? m ( S ))W m ( S ) + e (7Z - X )( s ^( s »7l'( S ,i? m ( S ))E m ( S ) 



c?s. 



Inserting the equations for A, //, A and // yields 

w m (t)-W m (t) = 



L 



(/i-A)( S ,-R m (s)) r (p-\)(s,R m (s)) t 



R 3 m(s)Em(s) R 3 m ( S )E m ( S ) 

R 2 rn (s) m R U.s) 

M^ +X){s ' Rm(s)) (w m ( S ) 3 ( S ,R„M)~E m (s)p( S ,R m ( S )) 



ds 



AttR, 



-47ri? m ( s ) e ^ +A )( s ^-( s »(iy m ( s )Xs,i? m ( s ))-£; m ( s )p( s ,i? m ( s )) 



: -Ft + -F 2 + F3 ; 



we have changed the order of terms to group the un-integrated source terms 
together. Using the same calculations as in the lemmas above we find 



|f x | < cj {\\ti*)-rtB)\\oo + \\H*)-H')\\°o 

+\\R(s)-R(s)\\ + \\W(s) -W(s)\\ ) ds, 



|F 2 | < CJ (^||M(s)-7J(s)||o + ||A(s)-A(s)|| 00 + ||m(s)-m(s)|| 00 

+ \\R(s)-R(s)\\ + \\W(s)-W(s)\\)ds, 
i? m ( S )e^ +A )( s ^-( s »(Ty m ( S )j( S ,i? m ( S ))--B m ( S )p( S ,i? m ( S ))) 



|F 3 | < An 



+ 47T 



R m (s)e {t * +X){S ' R ^ S » (w m (s)j(s,R m (s)) - E m (s)p(s,R m (s)) 



W m (s)(j(s,R m (s))-j(s,Rm(s))) 



R m (s)e^ +x){Ms)) 



-E m (s) (p(s,R m (s))-p(s,R m (s))) 



ds 



ds 
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-47T 



W m ( S )(j(s,Rm(s))-j(s,Rm(s))) 

-E m (s) (p(s,R m (s))-p(s,R m {s))) 



ds 



Denote the last term — the one with the differences j — j and p — p — by F4. Then 



\F«\ <C 



Ms) -Jl(s) Hoc + 1| A(s) - A(«)||oo + \\R(s) - R(s) 



+ \\W(s) - W(s)\\ + \\p(s) -p(a)||oo + Us) -=j(s) 



da+|F 4 |. 



We insert the definitions of the source terms into define k m (s):— 

e 2 H s > R rn(s)) E m [s)/ R m (s), recall the definitions of g m , n etc. from the proof 



of Lemma 3.13 and note that W m9m/ 'E m = Rm- Then 
SJo 



s Jo±^ E m (s) V 



< 



n=l 



1 ' "' is)f W n M n H^ - W n M n H^- ) (s) 



W — x W 2 x , 



~ / ^fc m ( S )M„( S )ij4 in ( S 

J / fern (s) =^J-TVw( s )H^ M (s)(W n M n - W n ~M n ) ( S )ds 

oJof-^ E m (s) 



n=l 



1 , W W Z \ 

g J o Y. k ^ s )9^) H kn^)[^M n -^M n \{s)ds 
1 r t N W 2 ( \ 

- / Y, k rn^)=^M n ( S )(g w -g TT )(s)(H^-H^ n )(s)d S 

1 r* N — — 

7 / y2k n (s)W n {s)M n {s) 

6 J n=l 

(H± n (s) (R m ~R n ) (s)-H^(s) (%n-X) (s))ds 



=: F 5 + F 6 + F 7 + F s + F 9 . 

Let Ii(s):={ne/fl\R n (s)-R m (s)\< 6} , I 2 (s) := {neM\\R n (s) -R m (s)\<6}. 
Then by Lemma 3J3, |ii(s)| < Ce~ 3 S and Xmeia(s) -^«( s ) < C(5. Hence, 



\F* 



\\<C ([[M*)-P(*)ll~ + ||A(a)-A( f 
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+ \\R(s)-R(s)\\ + \\W(s)-W(s)\\)ds, 
\ F e\<Cjf (\W„(s)-W n (s)\+e- 3 \M n (s)-M n (s)\\ds 

< cj^ (\\W(s)-W(s)\\ + \\M(s)-M(s)\\)ds, 
\F 7 \<Cj( Yl (\Rn-Rn\ + \W n ~W n \+e- 3 \M n ~M n \y s )ds 

< C [\\\R( S )-R(s)\\ + \\W( S )-W( S )\\ + \\M(s)-M(s)\\)ds, 



^-jj ^ M„{s)(\R m -R n \mm{\R n -R n \/S,l}Ys)ds 



n£/i(s)U/ 2 (s) 



< 



c 



M n (s)((5+\\R- R\\)mm{\\R- R\\/5,l}y s )ds 



n6/i(s)U/ 2 (s) 



< C / || R{s)-R{s) || ds. 
Jo 



Finally, let g n , m {s) :— k m (s)W n (s)M n (s). By Lemma \M n (s)\ <Ce 3 , hence 

||5n,m||oo<C£ 3 , 1 1 t/„, m || oo < Ce 3 , n,meJV. 
Thus Lemma [3.12 yields 



\F 9 \ < Ce- 3 max (|| fffciI |U + \\g k ,i\\oo) \\R(t) - R(t)\\ < C\\R(t) - R(t)\\. 



Collecting all these estimates and recalling Lemma 3.7 completes the proof. □ 

We are finally ready to prove our main results: 
Proof of Theorem 2.2 and Theorem 



3.2 



By Lemma 3.13 



\\M(t)-M(t)\\ < C[\\R{t)-R(t)\\+A 



where 



+ C j [\\W(s)-W(s)\\ + \\M(s)-M(s)\\)ds. 



A:= sup (||p(t)-?(*)||oo + lb(*)-P(t)||oo + ||j(*)-X*)llc 
te[o,r e , s ] v 



Combining this with Lemma 3.10 yields 



\\M(t)-M{t)\\ < C[\\R(t)-R{t)\\+A + e 



+ C J [\\W(s)-W(s)\\ + \\M(s)-M(s)\\jds. 
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Next, by Lemma 3.14, 



\\W(t)-W(t)\\ <c(\\R(t)-R(t)\\ + A) 

+ [\\\W(s)-W(s)\\ + \\M(s)-M(s)\\) 



ds. 



Adding the last two estimates and applying a Gronwall argument implies 



\W(t) - W(t)\\ + \\M(t) - M(t)\\ < Cl \\R(t) - R(t)\\ + A + e 



(3.5) 



Using Lemma 3.11 and Lemma 3.7 we find 



\Rn(t)-R n {t)\ <C [\Rn{s)-Rn(s)\ + \W n (s)-W n (s)\ 



+ \\n{s) -7i(s)||oo + 1| A(s) - A(s)||oo ds 



< CA 



-C 



(\\R(s)-R(s)\\ + \\W(s) - W(s)\\ + \\M(s)-M(s)\\) 



ds 



for all t E[0,T e .s] and n£j\f. Taking the maximum over n£j\f and inserting 
O) yields 



\\R(t)-R(t)\\<C[A + e + J \\R(s)-R{s)\\ds 

which via Gronwall implies that 

\\R(t)-R(t)\\<C{A + e). 



Now we apply Lemma 3.8 
A<C 



S + e/S in case of Theorem 3.2, 
5 2 + e/S in case of Theorem 2.2. 



(3.6) 



(3.7) 



Inserting this into (3.E) proves the corresponding error estimates on \\R(t)~ 
R(t)\\ in the two theorems^ The error estimates for || W(t) — W(t) \\ and 
\\M{t)-M{t)\\ follow from Q, and those for \\/j,(t) -7t(t)||oo, ||A(t) -A(t)||oo, 
and ||m (t) —mj^Wcg follow from Lemma 3.7. Under the assumption of Theo- 
rem 2fl the estimates for \\p(t) —'p(t)\\ 00 , \\p(t) — p(t)||oo and \\j(t) — X*)l|oo are 
valid due to Lemma |3.8j and Lemma 3.£. 

It remains to establish the assertions on the length T £ j of our approximation 
interval. Let C denote a constant for which the error estimates in Theorem L'.2| 



let S be so small that 



or Theorem |3.2| hold; note that C is independent of S and e. For Theorem 3.2 

e 1 
+ 5- 32D 2 C 
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Then 



Rn (t) < Rn (t) + - || < D + C (<J + |) < 

1 - 1 ^ R ^-J^<D + 2D^c(s + i)< 3 -D, 



R n (t) R n {t) R n (t)R n (t) ' ' V SJ-2 
\W n (t)\ < \W n {t)\ + \\W(t)-W(t)\\<D + c(6+^)<~D, 



\M n {t)\ < \M n (t)\+e 3 \\M(t)-M(t)\\<e 3 (D + c(t+^ 

e 2X(t,r) = e 2X(t,r) + £ ( _ _ s ( A+ A) (t ,r) < ^ + 16 ^ C fj + ^ < ^ 

r V (5/ 2 



for all t G [0,T E)5 ] and r > since m(f,r) = m(t,r) = for r < 1/(4D). Thus T ej<5 = 
T since otherwise the approximation interval could be extended beyond T e g. 



For Theorem 2.2. let S be so small that 



, e 1 e D 

5+ ^32W A5+ ^"2C- 

Then the above inequalities hold for all tG [0,r ei a], and in addition 

l|pWI|oo<||p(t)||oo + ||p(*)-p(*)[[ o<i3 + C , (d + ^)<|l). 

With the same argument as above, T Et g = T. This completes the proof of the 
main results. 



4 The fully discretized approximation 

In setting up a discretized version of a system like the Vlasov-Einstein system 
discretizing the phase space is the major step. For example, in several of the 
papers on numerical schemes for the Vlasov-Poisson or Vlasov-Maxwell systems 
mentioned in the introduction only this step is analyzed. However, due to the 
presence of the un-integrated source terms in the characteristic system, we can 
not discretize our system in time as easily . We have to modify the evolution 
equations slightly in order to get a convergence result, the idea being to discretize 



after integrating by parts in Lemma 3.12. Recall that the latter maneuver was 



essential for the convergence proof of the semi-discretized approximation, and 
we have to set up things in such a way that there is a substitute for this in the 
fully discretized case. 



Throughout this section the assumptions of Theorem |2.2| are to be satis- 
fied. Moreover, we assume e<8 2 and d<l/(8D). Let 0<r^;l be a small 
parameter, the time step. Before we discretize the evolution equations for 
Rn(t), W n (t), M n (t) we regroup the terms in these equations as follows; only 
the equations for W n (t) and M n (t) need to be considered, the discretization of 
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the equation for R n {t) being trivial: 



W n {t) 



(JZ-X)(t,R n (t)) _ 



Rl(t)E n (t) 



_ e (F+A)(i,K n (t)) m ( t -> R n{t)) 
R 2 n(t) 



E n (t) 



e 2X(t,R n (t)) _ 2 



f _ - 4tt< (i)e<"- A ^*" (*» (W„ (t)j(t, i?„ (t)) - £?„(*)?>(*, ^ (*))) 
H n {t) 



= e (J[-\)(t,Rn(t)) . 



L n c (7T+\)(t,R„(t)) m (^ R n('t)) ^ 

Tt n {t)E n {t) Tfjt) 

2X(t,R n (t)) 1 _ _ d - - 

+ - Mm (*) TA Rn w - Rm 



+■ 



e 2A(t,fl„(t)) x W 



VM m (t)T^l(i)|4rx.(i?nW-i? m W) 

" m= i Em(t) 

^(Jl-X)(t,R m (t)) _ e (Jl-X)(t,R n (t))^ 



--: FW n (t) + FW n (t) + FW n (t); 



note that 

j t X(Rn{t)-Rm(t)) = - 5 Xs(Rn(t)-R m {t)) 



e (p-\)(t,R n (t)) Wn(t) _ e (ji-X)(t,R m {t)) W m (t) 



E n {t) 



E m {t) 



Similarly, 
~M n {t) = M n {t) 

„2A(t,ft„(t)) 



(7T+ A) (t,R n (t)) m(t,R n (t)) W n {t) 

7f n (t) E n (t) 



+ ■ 



Rn(t) 



■ 2 '- (7T-A)(t,fi„(t)) 



-47ri?„(i)e 



J(t,i2„(t))-^^p(t,i2n(t)) 



M„(t) 



.,fil+AHt,fl»fO) m ( f ' fl "( f )) jM*) 



i£(t) £n(t) 

2X(t,I„(t)) i _ _ j _ _ 

■-^M m (t)£ m (t)- X (]? n (t)-iJ m (t)) 



+- 



e 2X(t,R n (t)) j w 



77 (/) ^^w^ww^w-^w) 

n V / m=l 







e (7Z+A)(t,fl n (t)) _ e (7Z+A)(t,K m (t)) N 



= : M„(t) 



FM n (t) + FM n (t)+™ n (t) 
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The terms FW n (t), FW n (t), FW n (t) and FM n (t), FMji), FM n (t) have 
natural discretizations given below; note that A\(R n (t) — R m (t)) is much bet- 
ter to discretizc and analyze than the term ^ (^R n — R„^j Xs(Rn{t) — R m (t)). 
Discretizing the equations above leads to the following Euler-like scheme: 

R —R .+ re (^-^)(«-.»)^i 

_j e 2X i (R n , i )_ 
Rn 



W n , i+1 := W n , l+ r FW n , t + ^= FW n<i + FW n<i , 



M„, <+ i := M„,, + tM„, 4 FM„, t + ^= FM ni + FM ni , 

cw 1 ._ e (^-A,)(fl„, I ) hn Ju,+X,)(R n ,) m ^ R ^) -p . 

r yv n,i ■— e 3 e 2 Ctt ^i 

FWl ■= 4^ E ^,«^n,i ~ ^ 

Rn,z m=1 T 

, e 2A((R n , 4 ) 1 2 1? . _ 

fW n>i := ^M^^^x^^-i?^) 

^n,i rn — 1 ^m.i 

AJii-XiKR^i) _ e (F i -Ai)(H„,i)^ i 

— —2 e 2A,(_R„,i) x(-Rn,i+l --Rm,»+l) ~ X(-RrM ~ Rm,i) 
tM n,i ■= = 2^ M m,i^m,i , 



Rn.i 



m—1 
N 



3 e 2A((R n , 4 ) j _ _ _ 

FM n i := -=2^M rn>< W rni< Xi(i2 n ,i-i2 m ,i) 

^n.* m=l 

^(^-AOCHn.i) _ e (7Ii-Ai)(iim,i)^ i 

with initial data 

(S n ,o,W n , ,Mn,o) := (fl n (0),F n (0),M n (0)) ) nGjV. 
We have abbreviated 



note that i?„,i+i can be computed without knowing the terms FW n i and FM n i 
where R n ,i+i appears in the difference quotient. The definitions of the approx- 
imating source terms etc. at the z-th time step are completely analogous to the 



34 



definitions of these quantities at time t = it in case of the phase space discretiza- 
tion. For example, for r > 0, 

1 r — - - 

Pi ( r ) •= M n,i E n,iX5 (r ~ Rn,i), 

n—1 

and analogously for p i and fa. Inserting these into the formulas for m(t), A(t), 
and Jl(t) defines corresponding approximations at the i-th time step of these 
quantities. 

To measure the errors for the fully discretized system we define 
ef := max \R n {tk) - Rn,k\, 
ef := max \W n (t k )-W n J, 
ef:= max e" 3 \~M n {t k ) - I„, t I ; 

note that these are monotonous in i. The convergence result for the fully dis- 
cretized scheme is as follows: 



Theorem 4.1 Suppose the assumptions of Theorem 2.& hold, e<S 2 , and let S 
and t be sufficiently small. Then 

ef + ef + ef < Ct, 
\\m(ir)-m l \\ 00 + \\jl(iT)-JI t \\ 0Q + \\X(iT)~Xi\\ 0C < Cr, 

IIp(*t) — Pilloo + lb(*-r) — Pilloo + ||j(i-r) — JiHoo < c- 

for all z G IN such that it 6 [0, T] , the time interval on which the solution f is to 
be approximated. 

As to the proof only some indications are given since once the fully dis- 
cretized system is set up in the appropriate way, what remains is to establish a 
series of auxiliary results most of which have analogues in the semi-discretized 
approximation which were proven in that situation. The detailed proof can be 
found in Let 



{ieJN|irG[0,T M ]}. 



The approximations are calculated as long as i£j e ^s and for all n£jV and 
k € {0, . . . ,i} the estimates 

£~ 3 M n , fc , |W n , fc |,i?„, fc , -L He 55 * Hoc, Halloo <4£ 

Rn,k 



hold; observe the analogy to (p.l|), fl3.3|). The set of all these indices i is denoted 
by J s .s,t For these i the error estimates from Theorem 4.1 are established. This 
is then used to show that if i £ Je,S,r and i + 1 £ J Et s then i + Jc,s,t, provided 
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the time step r is sufficiently small. Since T e .5 = T for sufficiently small 5 this 
completes this indication of the proof. 
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